Dynamical Decoupling Using Slow Pulses: Efficient Suppression of 1/f Noise 
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The application of dynamical decoupling pulses to a single qubit interacting with a linear harmònic 
oscillator bath with 1// spectral density is studied, and compared to the Ohmic case. Decoupling 
pulses that are slower than the fastest bath time-scale are shown to drastically reduce the decoher- 
ence rate in the 1/f case. Contrary to conclusions drawn from previous studies, this shows that 
dynamical decoupling pulses do not always have to be ultra-fast. Our results explain a recent experi- 
ment in which dephasing due to 1// charge noise affecting a charge qubit in a small superconducting 
elèctrode was successfully suppressed using spin-echo-type gate-voltage pulses. 
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The most serious problem in the physical implementa- 
tion of quantum information processing is that of main- 
taining quantum coherence. Decoherence due to inter- 
action with the environment can spoil the advantage of 
quantum algorithms 0] . One of the proposed remedies is 
the method of "dynamical decoupling" , or "bang-bang" 
(BB) pulses, in which strong and sufhciently fast pulses 
are applied to the system. In this manner one can ei- 
ther eliminate or symmctrize the system-bath Hamilto- 
nian so that system and bath are effectively decoupled 
ü H H \Ml U H H ES E3 The BB method was pro- 
posed in [3| , where a quantitative analysis was first per- 
formed for pure dephasing in the linear spin-boson model: 
Hsb = 9&z®B, where o~ z is the Pauli-z matrix and B is a 
Hermitian boson operator. The analysis was recently ex- 
tended to the non-linear spin-boson model, with similar 
conclusions about performance [3j. Decoupling also has 
been applied to the suppression of spontaneous emission 
0| and magnètic state decoherence induced by collisions 
in a vapor pj . Since the decoupling pulses are strong one 
ignores the evolution under Hs b while the pulses are on, 
and since the pulses are fast one ignores the evolution of 
the bath under its free Hamiltonian Hb during the pulse 
cycle. The latter assumption is usually stated as: 



Aí < 1/A 



uv i 



(1) 



where Aí is the pulse interval length and Ajjv is the 
high-frequency cutoff of the bath spectral density I(ui) 
[2| [see Eq. (J2J) below]. It can be shown that the overall 
system-bath coupling strength g is then renormalized by 
a factor AtAw after a cycle of decoupling pulses , and 
that the bath-induced error rate is reduced by a factor 
proportional to (AíA^y) 2 Q. A temperature T > sets 
an additional, thermal decoherence time scale that must 
be beat in order for the decoupling method to work 0, 0] • 
The conclusion (JIJ is extremely stringent, as the 
timescale Aí may be too small to be practically attain- 
able. Moreover, as we show below, and has been ar- 
gued before on the basis of the inverse quantum Zeno 
effect decoherence may be enhanced, rather than 



suppressed, if is not satisfied. Eq. is based on 
studies in which the bath was modeled as a system of 
harmònic oscillators, with a spectral density of the form 
I(lü) oc oj u e~ w / Kuv , with v > 0Lor using a fiat spectral 
density with a finitc cutoff Ajjv ]Éi , or without reference 
to a specific spectral density but emphasizing features of 
its high-frequency components 0, Q . However, a ubiqui- 
tous class of baths does not fall into this category, and we 
show here that then the condition Q is overly restrictive. 
This is the case for so-called 1/f noise, or more generally 
l/f a (a > 0). In these cases the bath spectral density 
decays as a power law, bounded between infrared (IR, 
lower) and ultraviolet (UV, upper) cutoffs Air and Ajjv, 
respectively. In quantum computer implementations this 
is often attributable to (but certainly not limited to) 
charge fiuctuations in elèctrodes providing control volt- 
ages. The need for such elèctrodes is widespread in quan- 
tum computer proposals, e.g., trapped ions (where ob- 
served 1/f noise was reported in E3)i quantum dots [l3| . 
doped silicon [14|. electrons on helium |lj| , and supercon- 
ducting qubits 16]. In the latter case, in a recent exper- 
iment involving a charge qubit in a small superconduct- 
ing elèctrode (Cooper-pair box), a spin-echo-type version 
of BB was successfully used to suppress low-frequency 
energy-level fiuctuations (causing dephasing) due tol// 
charge noise [ïtïj |. Here we explain the origin of such a 
result and discuss its general applicability. 

On the time scale í > 1/ ' Auv-, the details of the system- 
bath interaction and infernal bath dynamics become im- 
portant. These details are captured by the bath spectral 
density Since for í/f noise most of this density is 

concentrated in the low, rather than the high-end of the 
frequency range, it turns out that in this case BB with 
slow pulses (Aí > 1/Ajjv) depends more sensitively on 
the lower than on the upper cutoff. In particular, we 
show that the suppression of dephasing is more effective 
when the noise originates in a bath with 1/f spectrum 
than in the Ohmic case [v = 1 in Eq. J2J], owing to the 
abundance of IR modes in the formcr. In the following 
we present the results of our analysis contrasting BB for 



2 



1/f and Ohmic baths. 

Decoupling for spin-boson model. — We consider the 
linear spin-boson model including periòdic decoupling 
pulses. We first briefly review and somewhat simplify 
the results derived in [2|. We use kp = H = 1 units. The 
Hamiltonian is 



H = H. 



s 
e 



Hb + Hsb + Hp 

-f ^ u k b\b k + ^ o- z (g* k b k + 9kb\) + H f 



where the first (second) term governs the free system 
(bath) evolution; the third term is the (linear) system- 
bath interaction in which b k is the fcth-mode boson anni- 
hilation operator and g k is a coupling constant; and the 
last term is the fully controllable Hamiltonian generating 
the decoupling pulses: 



H P {t) = J2 V n (t) 



e ieto- í /2 (T e -ieta z /2 



n=l 



where the pulse amplitude V n (i) = V for t n < t < t n + t 
and otherwise, lasting for a duration r <C Aí, with t n = 
nAt being the time at which the nth pulse is applied. 
The properties of the bath are captured by its spectral 
density 



I ( UJ ) = ^2ò(uj -uj k )\g k \ 



(2) 



The reduced system density matrix is obtained from 
the total density matrix by tracing over the bath degrees 
of freedom 

Ps(t) = Tr fl [p(t)] = Tip [U(t)p s (0) <g> PB (0)U\t)] , 

where we have assumed a factorized initial condition be- 
tween the system and thermal bath, and U(t) is the time 

evolution generated by H: U(t) = Texp —i L dsH(s) 

(T denotes time ordering). We are interested in how 
decoupling improves the system coherence, defined as 
Poi(í) = (0| ps(t) In the interaction picture with re- 
spect to Hs and Hb the result in the absence of decou- 
pling pulses (free evolution) is: Pqi(*) = e_ro ^Voi(0)i 
where 



r (í) 



diu coth 



1 — COS LÜÍ2N 



m (3) 



(3 = l/(ItBT). In the Schròdinger picture there are 
oscillations at the natural frequency e, i.e., po\(t) — 



iet „I 



Similarly in the presence of the decoupling pulses, 
at t 2N = 2NAt, p^ihN) = e^e-W'l/jjiíO), 
where we can show from Eqs. (46), (47) of |2j that 



r P (A^,Aí) = 

4 / du coth 
The tan 2 



(4) 

/3lü\ l-COSUjt 2 N T , x , /wAí 

t) — ^ — Mtan v 

Aí) term (which was not found in 0) is 
the suppression factor arising from the decoupling proce- 
dure. In effect, the bath spectral density in the presence 
of decoupling pulses has been transformed from I(lü) to 
the singular spectral density I'(uj) — I{uj) tan 2 • 
Note, however, that the singularity of tan 2 at ojAt = 
(2n + í)tt for an integer n is canceled by the vanishing 
of 1 — coswÍ2Af at the same points, so Tp remains finite. 
Nevertheless, and as already pointed out in [2|, the value 
uiAt = 7r is special: In the limit N ^> 1 the integrand 
of Eq. Q is highly oscillatory for cuAt > ir, and grows 
to 16A^ 2 at u>At = 7r. Thus, decoherence suppression is 
effective when 



A uv At < 7T. 



(5) 



This is an upper bound on Aí that is independent of 
the specific form of I(uj). Note further that decoupling 
enhances decoherence from all modes with (4n+ l)7r/2 < 
ojAt < (4n + 3)7r/2, since for these vàlues tan 2 (wAí/2) > 
1. However, this effect may be quenched if the weight 
of these modes is sufficiently low; this is indeed what 
happens in the 1/f case. 

Results for 1/f and Ohmic spectral densities. — Let us 
now assume that the spectral density has the following 
form: 



I(ui) = jlü 1 ' 



±1, 



(6) 



with UV cutoff Kjjv and IR cutoff Ajp. Thus we are 
comparing 1/f noise (the case v = —1) to an Ohmic 
bath (the case v = 1, considered in 0]). 

To explain the effect of pulses qualitatively, we approx- 
imate tan 2 x by x 2 (l — 2x/ir)~ 1 , which allows us to obtain 
an explicit form for Tp for < KjjvAt < -k/2. We fur- 
ther expand coth ce w 1 + 2exp(— 2x) (x > 1). Then, the 
contribution to Tp for 1/f noise at low temperature is 
the sum of the zero temperature part 
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-,(T=0) 



(N,At) = 7 (Ai)' 



log 



A 



r/r 



log 



7T - A UV At 
7T - A/^Aí 



Ci (A ay í 2A r) + Ci (A t 2 jv) + O(Aí) 



(7) 



and the low temperature correction 



lf>°>(JV,At) = iM- 



log(l + T 2 í 2 



2AíT 



1 - 



1 + TH% N 



0(T 2 ) 



(8) 
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FIG. 1: Temporal behavior of the logarithm of the decoher- 
ence factors at T — 0. The initial coherence Poi(0) = 1. Pa- 
rameters are: 7 = 0.05, Auv = 10 for Ohmic and 7 = 0.25, 
Auv = 80 for 1//, Air = l,At = 0.025 for both. Thick sòlid 
(dashed) line: 1/f case with (without) decoupling pulses. 
Thin sòlid (dashed) line: Ohmic case with (without) decou- 
pling pulses. Eq. © was used for the case without decoupling 
pulses, while Eq. (J1J was used for the case with decoupling 
pulses at each t = tajv. The dotted line is our analytical result 
in Eq. 0. 

where Ci (Si) is the cosine (sine) integral. In Eq. JSJl, 
the limits Air — > and Ajjy — ¥ 00 are taken. All 
terms are finite in these limits. The first and second 
tornis in T^p ~~ ' ( N, At) (independent of Í2n) determine 

the asymptotic value T p (00, Aí); the remainder is a 
damped oscillatory part, given by the difference of two 
cosine integrals, that vanishes at long times. The sec- 
ond logarithmic term diverges as the pulse interval ap- 
proaches the inverse UV cutoff frequency time scale of 
the bath leading to decoherence enhancement from the 
tan 2 term in Eq. J3J. These behaviors are reflected in 
the exact solutions displayed in Fig-QJ The leading order 
finite temperature correction r^T^ (JV, Aí) can be sep- 
arated into two terms. The first term characterizes the 
asymptotic power law decay and the second term gives 
the initial damping and the asymptotic relaxation to the 
Í2jv-independent constant. 

In Fig.^the logarithm of the decoherence factors To(t) 
(free evolution) and T p (t) (pulsed evolution) for the 1/f 
and Ohmic cases are shown. The smaller is T, the more 
coherent is the evolution. The apparent oscillations with 



a frequency given by Auv are caused by the use of a 
sudden cutoff. Given the parameters used in Fig. 2] the 
Standard timescale condition Aí -c 1/Aj/y is not satis- 
fied in the 1/f case, while it is (AíAr/y = 0.25) in the 
Ohmic case. The most striking feature apparent in Fig.^ 
is the highly efficient suppression of decoherence in the 
case of 1/f noise, in spite of the seemingly unfavorable 
pulse interval length. In addition, it can be shown that 
decoherence due to the 1/f bath is accelerated when the 
IR cutoff is decreased, and is more sensitive to the IR 
cutoff than the Ohmic case. This is a direct consequence 
of the fact that most of the modes in 1/f spectrum are 
concentrated around Air. For 1/f baths we therefore 
cxpect slow and strong decoherence on a long time scale, 
that may be efficiently suppressed by relatively slow and 
strong pulses. A similar conclusion should be applica- 
ble to the more general class of baths with 1//" spectral 
density, since there too most of the bath spectral density 
is concentrated in the low frequency range. 

For our pure dephasing case at finite temperature, 
there is the thermal time scale tp = T -1 at which ther- 
mal fluctuations start affecting the system's coherence. 
In particular, for T 3> Ayy, decoherence is governed by 
the thermal fluctuations. In Fig. [21 a finite temperature 
result is shown. The decoupling pulses enhance the de- 
coherence for the Ohmic bath even at low temperatures, 
since for the parameters chosen the condition (1) is not 
satisfied. On the other hand, decoherence suppression in 
the 1/f case is highly effective. At high temperature, it 
has been argued on the basis of the Ohmic case, that de- 
coupling pulses faster than the thermal frequency T are 
required to suppress decoherence . Once again, the na- 
ture of the bath can qualitatively modify this conclusion. 
Thus decoupling by relatively slow pulses that obey the 
condition Aj/yAí ~ 1, can still be effective for decoher- 
ence suppression at elevated temperatures. However, as 
the temperature increases, the effective spectrum shifts 
toward low freqüències, and at the same time, the influ- 
ence of the environment increases. Overall, BB becomes 
ineffective irrespective of the type of bath. This explains 
the breakdown of decoherence suppression at T = 1000 
in Fig. |3| Note from the figure that the suppression of de- 
coherence for the 1/f bath is more effective than for the 
Ohmic bath throughout the whole temperature regime. 

For too slow pulses, BB accelerates the decoherence 
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|10| . For the Ohmic bath, as the interval approaches the 
threshold value Q from below, there is a crossover from 
decoherence suppression to decoherence enhancement, as 
shown in Fig. [3J For the í/f bath, suppression is still 
effective for longer pulse intervals as long as AtAuv < ^ 
is satisficd. 

It is of interest to compare our results with the gate 
voltage pulse experiment performed in [ïf| in a Cooper- 
pair box. The corresponding parameter vàlues in Eq. lf4"| 
are: 7 = 2Eq0 2 / 'e 2 h 2 , with the Josephson charging en- 
ergy E c = 122 \peV] and the constant a = (1.3 x 10~ 3 e) 2 
determined by the noise measurement. To achieve 90% 
decoherence suppression with A/fl = 100 [Hz] and Auv — 
10 [GHz] at k B T = 5 ]peV], the pulse interval Aí ~ 0.25 
[ns] is required from our analysis based on Eq. 10} with 
N = 1 . Although the pulse sequence of 0] differs from 
ours (theirs is the 7r/2 — tt — tt/2 spin-echo sequence), they 
play essentially the same role. Our Aí value roughly 
agrees with their value, Aí ~ 0.5 [ns], deduced from 
Fig. 2 in 0| . This agreement nicely illustrates the ex- 
perimental feasibility of BB in the case of 1// noise. The 
effectiveness of spin-echo type pulses in relation to su- 
perconducting qubits was also recently discussed in [l7j . 
The spin-boson model is appropriate for the study of 
1// noise due to a large number of weakly coupled back- 
ground charges[ÏÍ). 

Conclusions. — We have shown that the speed require- 
ment of the decoupling method should be stated rela- 
tive to the type of bath spectral density, and not just in 
terms of its upper cutoff (baths with bimodal distribu- 
tions provide another example of this 0). Most signifi- 
cantly, our exact results have demonstrated that BB can 
be expected to be highly effective in suppressing deco- 
herence due to the ubiquitous 1// noise, without having 
to satisfy the stringent time constraints that may ren- 
der the method overly difRcult to implement in other in- 
stances. We expect this to have significant implications, 
e.g., for suppression of noise due to charge fluctuations 
in elèctrodes providing control voltages in quantum com- 
putation. Such a result has already been obtained exper- 
imentally in a Cooper-pair box experiment |16| . and is 
predicted to apply to trapped-ion quantum computation 
as well [ïlj . 
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